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Abstract 

For the first time it is shown that for thin metallic films thickness of which 
does not exceed thickness of skin - layer, the problem of description of surface 
plasma oscillations allows analytical solution under arbitrary ratio of length of 
electron mean free path and thickness of a film. The dependance of frequency 
of surface plasma oscillations on wave number is deduced. We consider a case 
of specular - diffusive boundary value problems. 

Key words: degenerate collisional plasma, surface plasma oscillations, 
thin metallic film. 

PACS numbers: 73.50.-h Electronic transport phenomena in thin films, 
73.50.Mx High-frequency effects; plasma effects, 73.61.-r Electrical properties 
of specific thin films. 

Introduction 

Electromagnetic properties of metal films has been being a sub- 
ject of great interest for a long time already PQ [Bj. Problem of 
surface plasma oscillations has been a problem of special interest in 
recent time [7] [17] . It is connected as with theoretical interest to 
this problem, and with numerous practical appendices as well. Thus 
the majority of researches is founded on the description of properties 
of films with use of methods of macroscopical electrodynamics. Such 
approach is inadequate for thin films, since macroscopical electrody- 
namics is inapplicable for the description of films in the thickness of an 
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order of length of mean free path of electrons and less than this length. 
The electrons scattering on a surface demands kinetic consideration. 
It complicates the problem significantly. 

In the present work it is shown that for thin films, a thickness 
of which does not exceed a thickness of a skin - layer, the problem 
of description of surface plasma oscillations allows analytical solution 
under arbitrary ratio between length of mean free path of electrons 
and thickness of a film. The given work is a continuation of our work 
[18] in which the case when values of a magnetic field on top and 
bottom surface of films coincide was considered . Now the situation, 
when the signs of these values differ is investigated. 

Let us note, that the most part of reasonings carryied out be- 
low is true for more general case of conductive (in particular, semi- 
conductor) film. 

Problem statement 

Let us consider a thin metal film. 

We take Cartesian coordinate system with origin of coordinates on 
one of the surfaces of a slab, with axis x. directed deep into the slab 
and perpendicularly to the surface of the film. The axis z we will direct 
along the direction of propagation of the surface electromagnetic wave. 
We will note, that in this case magnetic field is directed along the axis 

y- 

Under such choice of system of coordinates the electric field vector 
and magnetic field vector have the following structure 

E = {E x (x 1 z,t),0,E z (x,z,t)}, H = {0,H y (x,z,t),0}. 

The origin of coordinates we will place at the bottom plane limiting 
a film. Let us designate a thickness of the film through d. 

Out of the film the electromagnetic field is described by the equa- 
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tions 

1 d E A T~l „ 

^" AE = 

and 

ia2H -AH = 0. 



cr dt 2 



,2 

Here c is the velocity of light, A is the Laplace operator. 
The solution of these equations decreasing at infinity point has the 
following form 

E = < . ' X<U ' (la) 

[ E 2e -^t+a(d-x)+zkz ^ X> d, 

and 

H = < ' X<U ' (16) 

[ H 2 e-^ +a ^^ +2fcz , x > d. 

Here cj is the wave frequency, k is the wave number, damping pa- 
rameter a is connected with these quantities by relation 



Ej and Hj (j = 1,2) are constant amplitudies. 

Further we search components of intensity vectors of electric and 
magnetic fields in the following form 

E x (x, z, t) = E x (x)e- M \ E z (x, z, t) = E z (x)e~ Mz , 

and 

Hy{x,z,t) = H y {x)e- lujt+lkz . 

Then the behaviour of electric and magnetic fields of the wave 
within the film is described by the following system of the differential 
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equations ([I] and [15] ) 



ifc-Ea, H # v = 0, 



dx c 

{ —E x - ikHy = —j x , (3) 
c c 

<ii7 w za; „ 47T 

-7^ + — ^ = — Jz- 

dx c c 

Here j is the current density. 

The equations (3) are satisfied out of the slab under the condition 
j = as well. 

Then impedance at the bottom surface of the layer (film) is defined 
as follows 

z0, = fS< i = 1 ' 2 ' (4) 

Quantities (j = 1,2) correspond to an impedance on the bot- 
tom layer surface. At the same time the quantity corresponds to 
magnetic field-symmetrical configuration of an external field. This is 
the case j ' = 1, when 

H y (0) = H y (d), E x (0) = E x (d), E z (0) = -E z {d). (5) 



This case has been considered in 
The quntity Z^ corresponds to configuration of an external field 
antisymmetric by magnetic field. It is the case j = 2, when 

H y (0) = -H y {d), E x {0) = -E x (d), £7,(0) = E z (d). (6) 

It is required to find a spatial dispersion of the surface plasmon, 
i.e. to find dependence of frequency of oscillations of eigen mode of 
the system (3) on quantity of the wave vector u = co(k). 

Let us consider a case when the width of the layer d is less than 
the depth of the skin - layer 5. We will note, that the depth of a skin 
- layer depends essentially on frequency of radiation, monotonously 



5 



decreasing in the process of growth of the last quantity. The quantity 
5 possesses minimum value in so-called infra-red case [2D] 8q = — , 

Up 

where u p is the plasma frequency. 

For typical metals [2D] ^0 ~ 10~ 5 cm. 

Thus for the films which thickness d is less 8q, our assumption is 
true for any frequencies. 

Quantities H y and E z vary a little on distances smaller than the 
depth of a skin - layer. Therefore at performance of the given assump- 
tion (d < 5q) thess fields will vary inside of slab. 



Surface plasmon. Antisymmetric configuration of magnetic 

field 

Let us consider the case 2 when £^(0) = E z (d). We can assume 
that in this layer z - projection of the electric field E z is constant. 
Then magnetic field change on the width of a layer can be defined 
from the third equation of the system (1) 

d 

H y (d) - H y (0) = ~—E z d + —I Jz (x)dx. (7) 



Thus 

j z {x) = a(x)E z , 

where a(x) is the conductivity which in general case depends on co- 
ordinate x. 

We introduce conductivity averaged by the slab thickness 

d d 

Od = jr^ J j z (x)dx = i J a(x)dx. 

Now we can rewrite the relation (7) in the form 

H y (d) - H y (0) = - l -^E z d+^^E z . 
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Considering symmetry of a magnetic field, from here we have 

Then for impedance (4) we have 

Z(2) = (8) 

uodyl + i -J 

In the same way as well as in (9), we receive 

= —. (9) 

to 

The dispersive equation for a surface plasma wave can be derived 
from the expressions (2), (8) and (9) 



2c Vc 2 /c 2 — uj 2 . 

ud + 4iri(Jdd u 

From the dispersive equation (10) we find the required spatial dis- 
persion 



k(uj) = — 
c 



4c 2 

l + ; 4^ —■ (11) 



Let us assume, that boundary conditions are specular - diffusive, p 
is the specularity coefficient. Let the relation kd <C 1 be true. Then in 
a low-frequency case, when u — > 0, the quantity o d can be presented 
in the form [H] 



and 



w d 

ad= $M a °> w= r (12 ^ 



1 3 , , f / 1 h 1- e~ wt , . 



$(w) w 2w 2 J U 3 t 5 Jl-pe 



Here / is the mean free path of electrons, p is the coefficient of 
specular reflection (specularity coefficient), do = cj 2 t/(47t) is the static 
conductivity of volume pattern, r = l/vp is the electron-transit time, 
vp is the Fermi velocity. 
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For arbitrary frequencies thee expression (12) and (13) will hold 
true under the condition, that it is necessary to use the following 
expression / — > (vft)/(1 — icur) as quantity /, and instead of ctq it's 
necessary to use the expression ctq — > <Jq/(1 — iur). 

Let us reduce the formula (11) to the form convenient for numerical 
calculations. We introduce the dimensionless parameters e = — and 



UJp 



Q = — . Then we can transform the formula (19) to the following 



form 



where 



Up 



oo 

— wt 



ip( W ) = 1 - —(l-p) (-- -)- -dt. 

rv 7 w v ' J \t 3 t 5 J 1 -pe- wt 

l 

In the case, when electrons reflect under specular condition from 
the film surface (i.e. at p = 1), the formula (14) becomes simpler and 
looks like 

k(u) = -Jl + ir . 2A2 , 15 

c W d^yu — no 2 + no 2 ) 2 

d 

where do = - is the half of a film thickness. 
2 

In dimensionless parametres the formula (15) can be written as 
follows 



^^y+mS^Sw- (i6) 

Under £ = we receive from here that 



It is clear that under e = from the last formula it follows that 
Im k(Q) = 0, i. e. in collisionless plasma plasmon damping is absent. 

From the formula (16) it is visible, that there exists such critical 
frequency uq = uo(e,d), that under co < co p Im k(co) > Re k(u), i.e. 
in field of subcritical frequencies sufrace plasmons do not exist. 
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We adduce the table of critical frequencies, referring to plasma 
(Langmuir) frequencies, Qo(e,d) = °^ — -, under e = 10 _1 for the 

Up 

case of specular boundary conditions (p = 1). 



Table 1 (critical frequencies) 



Film thickness d, nm 




1 


2 


3 


4 


Critical frequency 




0.101 


0.100 


0.097 


0.092 


Film thickness d, nm 


5 


6 


7 


8 


9 


Critical frequency 


0.086 


0.078 


0.067 


0.051 


0.023 



Let us adduce graphics on Figs. 1 - 8 of dependencies of real and 
imaginary parts of the wave vector on the ratio of frequencies uj/uj p 
under various values of frequency of electron collisions v, thickness of 
a film d and coefficient of specular reflection p. We will consider the 
case of sodium films, i.e. we take co p = 6.5 • 10 15 sec _1 , vp = 8.52 • 10 7 
cm/sec. 

Depending on quantity of parametres v, d, p the quantities Re k and 

Im k can essentially differ. So, at v = 10~ 5 cj p , d = 10 nanometer and 

p = 1 (Fig. 1) the quantity Re k surpasses Im k on some orders. 

Re k 

From Fig. 1 it is visible, that if to enter quantity Z = - — -, then 

z(o.i, nr 5 , io, i) = 2.i • io 4 , z(o.5, kt 5 , 10, l) = 3.8 • 10 4 . 

Besides, the quantity Re k is always positive under all values of 
parameters oj,is,d,p, while the quantity Im can be negative in the 
field of superhigh frequencies as well. 

Let us stop on existence of surface plasma waves (see Figs. 2 and 

3). 

Depending on quantities e, d,p two critical frequencies ujq and uji, 
such, that the inequality Im k < Re k is true under uo < cu < ui, and 
the inequality Im k > Re k is true under < to < cu p or u\ < u < u p 
can exist. 

In the last case of surface plasma waves do not exist. Let us consider 
the case v = 10 _1 w p ,j9 = 0.1. Results calculations we will present in 
the following table. 
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Table 2 (critical frequencies) 



Film thickness, nm 


First critical 


Second critical 


d 


frequency, uq 


frequency, uj\ 


1 


0.168 


0.904 


2 


0.130 


0.924 


3 


0.116 


0.929 


4 


0.107 


0.932 


5 


0.098 


0.934 


6 


0.089 


0.935 


7 


0.077 


0.935 


8 


0.063 


0.936 


9 


0.041 


0.936 


10 


0.000 


0.937 



The behaviour of the real and imaginary parts of a wave vector in 
dependence on a thickness of a film is presented on Figs. 4 and 5. The 
real part has a sharp maximum nearby the plasma resonance to ~ cu p . 

Let us note, the more is thickness of a film the less are values 
of the real part under each value of frequency of oscillations of an 
electromagnetic field. The imaginary part has the same behaviour. 
But its values are less than values of real parts significantly. Unlike 
the real part in the field of the superhigh frequencies the values of the 
imaginary part become negative. 

On Figs. 6 and 7 dependences of the real and imaginary parts 
of wave vector on quantity of collision frequencies of electrons are 
presented. The less is the quantity v, the more are the values of the 
real part. For the imaginary part inverse relation takes place. Namely, 
the less is the quantity of collision frequencies, the less is the value of 
an imaginary part. It means, that at electron collisions frequency 
increase the attenuation of surface plasma waves becomes stronger (in 
the field of subcritical frequencies). 

On Fig. 8 dependence of an imaginary part of the wave vector on 
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quantity of coefficient of specular reflection in the field of subcritical 
frequencies is presented. Graphics show, that with the growth of coef- 
ficient of specular reflection the values of an imaginary part decrease. 
It means, that damping of plasma waves by that becomes the stronger, 
the less are quantities of coefficient of specular reflection. 

Conclusion 

In the present work the dispersion relation for surface plasmons is 
deduced. We consider the case of an antisymmetric configuration of x 
- component of the electric field and y - component of magnetic field, 
and symmetric z - component of electric field. We consider the case 
of specular - diffusive boundary value conditions. 
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Figure 1. Real and imaginary parts of wave number, film thickness d = 10 nm, collision 
electron frequency v = 10~ 5 w p 1/sec, specular reflection, p = 1. 
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Figure 2. Real and imaginary parts of wave number, film thickness d = 100 nm, collision 
electron frequency v = 0.1w p 1/sec, coefficient of specular reflection, p = 0.1. 
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Figure 3. Real and imaginary parts of wave number, film thickness d = 100 nm, collision 
electron frequency v = 0.02w p 1/sec, coefficient of specular reflection, p = 0.1. 




Figure 4. Real and imaginary parts of wave number, collision electron frequency 
v = 10~ 3 oj p 1/sec, p = 0.5, curves 1, 2, 3, 4 correspond to values of the film thickness 
d= 10,25,50,100 nm. 
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Figure 5. Imaginary part of wave number, v = 10~ 3 u;p 1/sec, p = 0.5, curves 1, 2, 3, 4 
correspond to values of the film thickness d = 10, 25, 50, 100 nm. Subcritical frequencies: 
< oj < uj p 1/sec. 
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Figure 6. Real part of wave number, curves 1, 2, 3 correspond to values of electron 
collision frequencies v = 10 _5 oj p , 5 • lO -2 ^, lO -1 ^ 1/sec. Film thickness d = 10 nm, 
coefficient of specular reflection p = 1 (case of specular reflection). 
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Figure 7. Imaginary part of wave number, curves 1, 2, 3 correspond to values of electron 
collision frequencies v = 10 _5 Wp, 5 • 10 -2 u;p, 10 _1 w p 1/sec. Film thickness d = 10 nm, 
coefficient of specular reflection p = 1 (case of specular reflection). 
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Figure 8. Imaginary part of the wave number Im k(u). Curves 1, 2, 3 correspond to 
values of coefficient of specular reflection p = 0, 0.5, 1. Electron collision frequencies is 
equal v = 10~ 3 w p 1/sec. Film thickness d = 1 nm. 



